A large-effect fitness trade-off across environments is explained by a single mutation affecting cold acclimation.

Identifying the genetic basis of local adaptation and fitness trade-offs across environments is a central goal of evolutionary biology. Cold acclimation is an adaptive plastic response for surviving seasonal freezing, and costs of acclimation may be a general mechanism for fitness trade-offs across environments in temperate zone species. Starting with locally adapted ecotypes of Arabidopsis thaliana from Italy and Sweden, we examined the fitness consequences of a naturally occurring functional polymorphism in CBF2. This gene encodes a transcription factor that is a major regulator of cold-acclimated freezing tolerance and resides within a locus responsible for a genetic trade-off for long-term mean fitness. We estimated the consequences of alternate genotypes of CBF2 on 5-y mean fitness and fitness components at the native field sites by comparing near-isogenic lines with alternate genotypes of CBF2 to their genetic background ecotypes. The effects of CBF2 were validated at the nucleotide level using gene-edited lines in the native genetic backgrounds grown in simulated parental environments. The foreign CBF2 genotype in the local genetic background reduced long-term mean fitness in Sweden by more than 10%, primarily via effects on survival. In Italy, fitness was reduced by more than 20%, primarily via effects on fecundity. At both sites, the effects were temporally variable and much stronger in some years. The gene-edited lines confirmed that CBF2 encodes the causal variant underlying this genetic trade-off. Additionally, we demonstrated a substantial fitness cost of cold acclimation, which has broad implications for potential maladaptive responses to climate change.

Can heterosis and inbreeding depression explain the maintenance of outcrossing in a cleistogamous perennial?

PREMISE: What maintains mixed mating is an evolutionary enigma. Cleistogamy-the production of both potentially outcrossing chasmogamous and obligately selfing cleistogamous flowers on the same individual plant-is an excellent system to study the costs of selfing. Inbreeding depression can prevent the evolution of greater selfing within populations, and heterosis in crosses between populations may further tip the balance in favor of outcrossing. Few empirical estimates of inbreeding depression and heterosis in the same system exist for cleistogamous species. METHODS: We investigate the potential costs of selfing by quantifying inbreeding depression and heterosis in three populations of the cleistogamous perennial Ruellia humilis Nutt (Acanthaceae). We performed three types of hand-pollinations-self, outcross-within, and outcross-between populations-and measured seed number, germination, total flower production, and estimated cumulative fitness for the resulting progeny in a greenhouse experiment. RESULTS: We found moderate inbreeding depression for cumulative fitness (<30%) in two populations, but outbreeding depression for crosses within a third population (-26%). For between-population crosses, there was weak to modest heterosis (11-47%) in two of the population combinations, but modest to strong outbreeding depression (-21 to -71%) in the other four combinations. CONCLUSIONS: Neither inbreeding depression nor heterosis was of sufficient magnitude to explain the continued production of chasmogamous flowers given the relative energetic advantage of cleistogamous flowers previously estimated for these populations. Outbreeding depression either within or between populations makes the maintenance of chasmogamous flowers even harder to explain. More information is needed on the genetic basis of cleistogamy to resolve this conundrum.

PREMISA: Lo que mantiene los sistemas de apareamiento mixto aún es un enigma. La cleistogamia, la producción de flores con potencial de cruzamiento casmógamas, y de flores cleistógamas obligadamente autofecundadas en la misma planta, es un excelente sistema para estudiar los costos de la autofecundación. La depresión endogámica puede prevenir la evolución hacia una mayor autofecundación dentro de las poblaciones, y la heterosis de los cruces entre poblaciones puede inclinar aún más la balanza a favor del cruzamiento. Existen pocas estimaciones empíricas de depresión endogámica y heterosis en el mismo sistema para especies cleistógamas. MÉTODOS: Investigamos los costos potenciales de la autofecundación cuantificando la depresión endogámica y la heterosis en tres poblaciones de la perenne cleistógama Ruellia humilis Nutt (Acanthaceae). Realizamos autopolinizaciones manuales, y cruces dentro y entre poblaciones, Medimos el número de semillas, la germinación, la producción total de flores y estimamos la acumulación de fitness para la progenie resultante en un experimento de invernadero. RESULTADOS: Encontramos depresión endogámica moderada para fitness acumulado (<30%) en dos poblaciones, pero depresión exogámica para cruces dentro de la tercera población (-26%). Entre cruces de población, hubo heterosis de débil a modesta (11-47%) en dos de las combinaciones de poblaciones, pero depresión exogámica moderada a fuerte (-21 a -71%) en las otras cuatro combinaciones. CONCLUSIONES: Ni la depresión endogámica, ni la heterosis fueron de suficiente magnitud para explicar la producción continua de flores casmógamas dada la ventaja energética relativa de las flores cleistógamas previamente estimadas para estas poblaciones. La depresión exogámica, ya sea dentro o entre poblaciones, hace que el mantenimiento de las flores casmógamas sea aún más difícil de explicar. Se necesita más información sobre la base genética de la cleistogamia para resolver este enigma.

Ecological genetics of local adaptation in Arabidopsis: An 8-year field experiment.

There is considerable evidence for local adaptation in nature, yet important questions remain regarding its genetic basis. How many loci are involved? What are their effect sizes? What is the relative importance of conditional neutrality versus genetic trade-offs? Here we address these questions in the self-pollinating, annual plant Arabidopsis thaliana. We used 400 recombinant inbred lines (RILs) derived from two locally adapted populations in Italy and Sweden, grew the RILs and parents at the parental locations, and mapped quantitative trait loci (QTL) for mean fitness (fruits/seedling planted). We previously published results from the first 3 years of the study, and here add five additional years, providing a unique opportunity to assess how temporal variation in selection might affect QTL detection and classification. We found 10 adaptive and one maladaptive QTL in Italy, and six adaptive and four maladaptive QTL in Sweden. The discovery of maladaptive QTL at both sites suggests that even locally adapted populations are not always at their genotypic optimum. Mean effect sizes for adaptive QTL, 0.97 and 0.55 fruits in Italy and Sweden, respectively, were large relative to the mean fitness of the RILs (approximately 8 fruits/seedling planted at both sites). Both genetic trade-offs (four cases) and conditional neutrality (seven cases) contribute to local adaptation in this system. The 8-year dataset provided greater power to detect QTL and to estimate their locations compared to our previous 3-year study, identifying one new genetic trade-off and resolving one genetic trade-off into two conditionally adaptive QTL.

Life-history trade-offs and the genetic basis of fitness in Arabidopsis thaliana.

Resources allocated to survival cannot be used to increase fecundity, but the extent to which this trade-off constrains adaptation depends on overall resource status. Adaptation to local environmental conditions may therefore entail the evolution of traits that increase the amount of resources available to individuals (their resource status or 'condition'). We examined the relative contribution of trade-offs and increased condition to adaptive evolution in a recombinant inbred line population of Arabidopsis thaliana planted at the native sites of the parental ecotypes in Italy and Sweden in 2 years. We estimated genetic correlations among fitness components based on genotypic means and explored their causes with QTL mapping. The local ecotype produced more seeds per fruit than did the non-local ecotype, reflected in stronger adaptive differentiation than was previously shown based on survival and fruit number only. Genetic correlations between survival and overall fecundity, and between number of fruits and number of seeds per fruit, were positive, and there was little evidence of a trade-off between seed size and number. Quantitative trait loci for these traits tended to map to the same regions of the genome and showed positive pleiotropic effects. The results indicate that adaptive differentiation between the two focal populations largely reflects the evolution of increased ability to acquire resources in the local environment, rather than shifts in the relative allocation to different life-history traits. Differentiation both in phenology and in tolerance to cold is likely to contribute to the advantage of the local genotype at the two sites.

Genetic and physiological mechanisms of freezing tolerance in locally adapted populations of a winter annual.

PREMISE: Despite myriad examples of local adaptation, the phenotypes and genetic variants underlying such adaptive differentiation are seldom known. Recent work on freezing tolerance and local adaptation in ecotypes of Arabidopsis thaliana from Italy and Sweden provides an essential foundation for uncovering the genotype-phenotype-fitness map for an adaptive response to a key environmental stress. METHODS: We examined the consequences of a naturally occurring loss-of-function (LOF) mutation in an Italian allele of the gene that encodes the transcription factor CBF2, which underlies a major freezing-tolerance locus. We used four lines with a Swedish genetic background, each containing a LOF CBF2 allele. Two lines had introgression segments containing the Italian CBF2 allele, and two contained deletions created using CRISPR-Cas9. We used a growth chamber experiment to quantify freezing tolerance and gene expression before and after cold acclimation. RESULTS: Freezing tolerance was lower in the Italian (11%) compared to the Swedish (72%) ecotype, and all four experimental CBF2 LOF lines had reduced freezing tolerance compared to the Swedish ecotype. Differential expression analyses identified 10 genes for which all CBF2 LOF lines, and the IT ecotype had similar patterns of reduced cold responsive expression compared to the SW ecotype. CONCLUSIONS: We identified 10 genes that are at least partially regulated by CBF2 that may contribute to the differences in cold-acclimated freezing tolerance between the Italian and Swedish ecotypes. These results provide novel insight into the molecular and physiological mechanisms connecting a naturally occurring sequence polymorphism to an adaptive response to freezing conditions.

Heterosis is common and inbreeding depression absent in natural populations of Arabidopsis thaliana.

The importance of genetic drift in shaping patterns of adaptive genetic variation in nature is poorly known. Genetic drift should drive partially recessive deleterious mutations to high frequency, and inter-population crosses may therefore exhibit heterosis (increased fitness relative to intra-population crosses). Low genetic diversity and greater genetic distance between populations should increase the magnitude of heterosis. Moreover, drift and selection should remove strongly deleterious recessive alleles from individual populations, resulting in reduced inbreeding depression. To estimate heterosis, we crossed 90 independent line pairs of Arabidopsis thaliana from 15 pairs of natural populations sampled across Fennoscandia and crossed an additional 41 line pairs from a subset of four of these populations to estimate inbreeding depression. We measured lifetime fitness of crosses relative to parents in a large outdoor common garden (8,448 plants in total) in central Sweden. To examine the effects of genetic diversity and genetic distance on heterosis, we genotyped parental lines for 869 SNPs. Overall, genetic variation within populations was low (median expected heterozygosity = 0.02), and genetic differentiation was high (median FST  = 0.82). Crosses between 10 of 15 population pairs exhibited significant heterosis, with magnitudes of heterosis as high as 117%. We found no significant inbreeding depression, suggesting that the observed heterosis is due to fixation of mildly deleterious alleles within populations. Widespread and substantial heterosis indicates an important role for drift in shaping genetic variation, but there was no significant relationship between fitness of crosses relative to parents and genetic diversity or genetic distance between populations.

Combining population genomics and fitness QTLs to identify the genetics of local adaptation in Arabidopsis thaliana.

Evidence for adaptation to different climates in the model species Arabidopsis thaliana is seen in reciprocal transplant experiments, but the genetic basis of this adaptation remains poorly understood. Field-based quantitative trait locus (QTL) studies provide direct but low-resolution evidence for the genetic basis of local adaptation. Using high-resolution population genomic approaches, we examine local adaptation along previously identified genetic trade-off (GT) and conditionally neutral (CN) QTLs for fitness between locally adapted Italian and Swedish A. thaliana populations [Ågren J, et al. (2013) Proc Natl Acad Sci USA 110:21077-21082]. We find that genomic regions enriched in high FST SNPs colocalize with GT QTL peaks. Many of these high FST regions also colocalize with regions enriched for SNPs significantly correlated to climate in Eurasia and evidence of recent selective sweeps in Sweden. Examining unfolded site frequency spectra across genes containing high FST SNPs suggests GTs may be due to more recent adaptation in Sweden than Italy. Finally, we collapse a list of thousands of genes spanning GT QTLs to 42 genes that likely underlie the observed GTs and explore potential biological processes driving these trade-offs, from protein phosphorylation, to seed dormancy and longevity. Our analyses link population genomic analyses and field-based QTL studies of local adaptation, and emphasize that GTs play an important role in the process of local adaptation.

Genetic basis of photosynthetic responses to cold in two locally adapted populations of Arabidopsis thaliana.

Local adaptation is common, but the traits and genes involved are often unknown. Physiological responses to cold probably contribute to local adaptation in wide-ranging species, but the genetic basis underlying natural variation in these traits has rarely been studied. Using a recombinant inbred (495 lines) mapping population from locally adapted populations of Arabidopsis thaliana from Sweden and Italy, we grew plants at low temperature and mapped quantitative trait loci (QTLs) for traits related to photosynthesis: maximal quantum efficiency (Fv/Fm), rapidly reversible photoprotection (NPQfast), and photoinhibition of PSII (NPQslow) using high-throughput, whole-plant measures of chlorophyll fluorescence. In response to cold, the Swedish line had greater values for all traits, and for every trait, large effect QTLs contributed to parental differences. We found one major QTL affecting all traits, as well as unique major QTLs for each trait. Six trait QTLs overlapped with previously published locally adaptive QTLs based on fitness measured in the native environments over 3 years. Our results demonstrate that photosynthetic responses to cold can vary dramatically within a species, and may predominantly be caused by a few QTLs of large effect. Some photosynthesis traits and QTLs probably contribute to local adaptation in this system.

Adaptive divergence in flowering time among natural populations of Arabidopsis thaliana: Estimates of selection and QTL mapping.

To identify the ecological and genetic mechanisms of local adaptation requires estimating selection on traits, identifying their genetic basis, and evaluating whether divergence in adaptive traits is due to conditional neutrality or genetic trade-offs. To this end, we conducted field experiments for three years using recombinant inbred lines (RILs) derived from two ecotypes of Arabidopsis thaliana (Italy, Sweden), and at each parental site examined selection on flowering time and mapped quantitative trait loci (QTL). There was strong selection for early flowering in Italy, but weak selection in Sweden. Eleven distinct flowering time QTL were detected, and for each the Italian genotype caused earlier flowering. Twenty-seven candidate genes were identified, two of which (FLC and VIN3) appear under major flowering time QTL in Italy. Seven of eight QTL in Italy with narrow credible intervals colocalized with previously reported fitness QTL, in comparison to three of four in Sweden. The results demonstrate that the magnitude of selection on flowering time differs strikingly between our study populations, that the genetic basis of flowering time variation is multigenic with some QTL of large effect, and suggest that divergence in flowering time between ecotypes is due mainly to conditional neutrality.

Genetics of water use physiology in locally adapted Arabidopsis thaliana.

Identifying the genetic basis of adaptation to climate has long been a goal in evolutionary biology and has applications in agriculture. Adaptation to drought represents one important aspect of local adaptation, and drought is the major factor limiting agricultural yield. We examined local adaptation between Sweden and Italy Arabidopsis thaliana ecotypes, which show contrasting levels of water availability in their local environments. To identify quantitative trait loci (QTL) controlling water use physiology traits and adaptive trait QTL (genomic regions where trait QTL and fitness QTL colocalize), we performed QTL mapping on 374F9 recombinant inbred lines in well-watered and terminal drought conditions. We found 72 QTL (32 in well-watered, 31 in drought, 9 for plasticity) across five water use physiology traits: δ(13)C, rosette area, dry rosette weight, leaf water content and percent leaf nitrogen. Some of these genomic regions colocalize with fitness QTL and with other physiology QTL in defined hotspots. In addition, we found evidence of both constitutive and inducible water use physiology QTL. Finally, we identified highly divergent candidate genes, in silico. Our results suggest that many genes with minor effects may influence adaptation through water use physiology and that pleiotropic water use physiology QTL have fitness consequences.

Factors influencing the effect size distribution of adaptive substitutions.

The distribution of effect sizes of adaptive substitutions has been central to evolutionary biology since the modern synthesis. Early theory proposed that because large-effect mutations have negative pleiotropic consequences, only small-effect mutations contribute to adaptation. More recent theory suggested instead that large-effect mutations could be favoured when populations are far from their adaptive peak. Here we suggest that the distributions of effect sizes are expected to differ among study systems, reflecting the wide variation in evolutionary forces and ecological conditions experienced in nature. These include selection, mutation, genetic drift, gene flow, and other factors such as the degree of pleiotropy, the distance to the phenotypic optimum, whether the optimum is stable or moving, and whether new mutation or standing genetic variation provides the source of adaptive alleles. Our goal is to review how these factors might affect the distribution of effect sizes and to identify new research directions. Until more theory and empirical work is available, we feel that it is premature to make broad generalizations about the effect size distribution of adaptive substitutions important in nature.

Increased heterosis in selfing populations of a perennial forb.

Quantifying the importance of random genetic drift in natural populations is central to understanding the potential limits to natural selection. One approach is to estimate the magnitude of heterosis, the increased fitness of progeny derived from crosses between populations relative to crosses within populations caused by the heterozygous masking of deleterious recessive or nearly recessive alleles that have been fixed by drift within populations. Self-fertilization is expected to reduce the effective population size by half relative to outcrossing, and population bottlenecks may be common during the transition to selfing. Therefore, chance fixation of deleterious alleles due to drift in selfing populations should increase heterosis between populations. Increased homozygosity due to fixation or loss of alleles should also decrease inbreeding depression within populations. Most populations of the perennial herb Arabidopsis lyrata ssp. lyrata are self-incompatible (SI), but several have evolved self-compatibility and are highly selfing. We quantified heterosis and inbreeding depression in two predominantly self-compatible (SC) and seven SI populations in a field common garden experiment within the species' native range and examined the correlation between these metrics to gauge the similarity in their genetic basis. We measured proportion germination in the lab, and survival and fecundity (flower and seed production) for 2 years in the field, and calculated estimates of cumulative fitness. We found 7.2-fold greater heterosis in SC compared with SI populations, despite substantial heterosis in SI populations (56 %). Inbreeding depression was >61 %, and not significantly different between SC and SI populations. There was no correlation between population estimates of heterosis and inbreeding depression, suggesting that they have somewhat different genetic bases. Combined with other sources of information, our results suggest a history of bottlenecks in all of these populations. The bottlenecks in SC populations may have been severe, but their strong inbreeding depression remains enigmatic.

QTL mapping of freezing tolerance: links to fitness and adaptive trade-offs.

Local adaptation, defined as higher fitness of local vs. nonlocal genotypes, is commonly identified in reciprocal transplant experiments. Reciprocally adapted populations display fitness trade-offs across environments, but little is known about the traits and genes underlying fitness trade-offs in reciprocally adapted populations. We investigated the genetic basis and adaptive significance of freezing tolerance using locally adapted populations of Arabidopsis thaliana from Italy and Sweden. Previous reciprocal transplant studies of these populations indicated that subfreezing temperature is a major selective agent in Sweden. We used quantitative trait locus (QTL) mapping to identify the contribution of freezing tolerance to previously demonstrated local adaptation and genetic trade-offs. First, we compared the genomic locations of freezing tolerance QTL to those for previously published QTL for survival in Sweden, and overall fitness in the field. Then, we estimated the contributions to survival and fitness across both field sites of genotypes at locally adaptive freezing tolerance QTL. In growth chamber studies, we found seven QTL for freezing tolerance, and the Swedish genotype increased freezing tolerance for five of these QTL. Three of these colocalized with locally adaptive survival QTL in Sweden and with trade-off QTL for overall fitness. Two freezing tolerance QTL contribute to genetic trade-offs across environments for both survival and overall fitness. A major regulator of freezing tolerance, CBF2, is implicated as a candidate gene for one of the trade-off freezing tolerance QTL. Our study provides some of the first evidence of a trait and gene that mediate a fitness trade-off in nature.

Flowering time QTL in natural populations of Arabidopsis thaliana and implications for their adaptive value.

The genetic basis of phenotypic traits is of great interest to evolutionary biologists, but their contribution to adaptation in nature is often unknown. To determine the genetic architecture of flowering time in ecologically relevant conditions, we used a recombinant inbred line population created from two locally adapted populations of Arabidopsis thaliana from Sweden and Italy. Using these RILs, we identified flowering time QTL in growth chambers that mimicked the natural temperature and photoperiod variation across the growing season in each native environment. We also compared the genomic locations of flowering time QTL to those of fitness (total fruit number) QTL from a previous three-year field study. Ten total flowering time QTL were found, and in all cases, the Italy genotype caused early flowering regardless of the conditions. Two QTL were consistent across chamber environments, and these had the largest effects on flowering time. Five of the fitness QTL colocalized with flowering time QTL found in the Italy conditions, and in each case, the local genotype was favoured. In contrast, just two flowering time QTL found in the Sweden conditions colocalized with fitness QTL and in only one case was the local genotype favoured. This implies that flowering time may be more important for adaptation in Italy than Sweden. Two candidate genes (FLC and VIN3) underlying the major flowering time QTL found in the current study are implicated in local adaptation.

Genetic mapping of adaptation reveals fitness tradeoffs in Arabidopsis thaliana.

Organisms inhabiting different environments are often locally adapted, and yet despite a considerable body of theory, the genetic basis of local adaptation is poorly understood. Unanswered questions include the number and effect sizes of adaptive loci, whether locally favored loci reduce fitness elsewhere (i.e., fitness tradeoffs), and whether a lack of genetic variation limits adaptation. To address these questions, we mapped quantitative trait loci (QTL) for total fitness in 398 recombinant inbred lines derived from a cross between locally adapted populations of the highly selfing plant Arabidopsis thaliana from Sweden and Italy and grown for 3 consecutive years at the parental sites (>40,000 plants monitored). We show that local adaptation is controlled by relatively few genomic regions of small to modest effect. A third of the 15 fitness QTL we detected showed evidence of tradeoffs, which contrasts with the minimal evidence for fitness tradeoffs found in previous studies. This difference may reflect the power of our multiyear study to distinguish conditionally neutral QTL from those that reflect fitness tradeoffs. In Sweden, but not in Italy, the local genotype underlying fitness QTL was often maladaptive, suggesting that adaptation there is constrained by a lack of adaptive genetic variation, attributable perhaps to genetic bottlenecks during postglacial colonization of Scandinavia or to recent changes in selection regime caused by climate change. Our results suggest that adaptation to markedly different environments can be achieved through changes in relatively few genomic regions, that fitness tradeoffs are common, and that lack of genetic variation can limit adaptation.

Small effective size limits performance in a novel environment.

Understanding what limits or facilitates species' responses to human-induced habitat change can provide insight for the control of invasive species and the conservation of small populations, as well as an arena for studying adaptation to realistic novel environments. Small effective size of ancestral populations could limit the establishment in, or response to, a novel or altered habitat because of low genetic variation for ecologically important traits, and/or because small populations harbor fixed deleterious mutations. I estimated the fitness of individuals from populations of the endangered plant Hypericum cumulicola, of known census and effective size, transplanted into native scrub habitat and unpaved roadsides, which are a novel habitat for this species. I found a significant positive relationship between estimates of population size and mean fitness, but only in the novel roadside habitat. Fitness was more than 200% greater in the roadside habitat than the scrub, mostly due to increased fecundity. These results combined with previous estimates of heterosis in this species suggest that fixed deleterious mutations could contribute to lower fitness of field transplants from small populations in the novel environment.

Effects of population size and isolation on heterosis, mean fitness, and inbreeding depression in a perennial plant.

• In small isolated populations, genetic drift is expected to increase chance fixation of partly recessive, mildly deleterious mutations, reducing mean fitness and inbreeding depression within populations and increasing heterosis in outcrosses between populations. • We estimated relative effective sizes and migration among populations and compared mean fitness, heterosis, and inbreeding depression for eight large and eight small populations of a perennial plant on the basis of fitness of progeny produced by hand pollinations within and between populations. • Migration was limited, and, consistent with expectations for drift, mean fitness was 68% lower in small populations; heterosis was significantly greater for small (mean = 70%, SE = 14) than for large populations (mean = 7%, SE = 27); and inbreeding depression was lower, although not significantly so, in small (mean = -0.29%, SE = 28) than in large (mean = 0.28%, SE = 23) populations. • Genetic drift promotes fixation of deleterious mutations in small populations, which could threaten their persistence. Limited migration will exacerbate drift, but data on migration and effective population sizes in natural populations are scarce. Theory incorporating realistic variation in population size and patterns of migration could better predict genetic threats to small population persistence.